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Abstract. We describe the use of array expressions as constraints, 
which represents a consequent generalisation of the element constraint. 
Constraint propagation for array constraints is studied theoretically, and 
for a set of domain reduction rules the local consistency they enforce, arc- 
consistency, is proved. An efficient algorithm is described that encapsu- 
lates the rule set and so inherits the capability to enforce arc-consistency 
from the rules. 



1 Introduction 

Many problems can be modelled advantageously using "look up" functionality: 
access an object given an index. Imperative programming languages provide ar- 
rays for this. With i one of 1, 2, 3 and a definition such as integer a [3], 
the construct a[i] represents an integer variable, while with a definition 
b[] = {5, 7, 9} the 'value' of i according to table b can be looked up by 
x = b[i]. 

A usual condition for look-up expressions to be valid is that the index be 
known when the expression is evaluated. In a constraint programming environ- 
ment this is a restriction that can be disposed of. The binary element constraint 
(originally in CHIP, Q), semantically equivalent to a lookup expression using a 
1-dimensional array, allows a variable as the index and a variable for the result, 
constraining both. It has proved very beneficial to allow a variable for the in- 
dex. Many important problems (scheduling, resource allocation, etc.) modelled 
as CSPs make use of this constraint. 

OPL, a modelling language for combinatorial optimisation and constraint 
programming (Jll|), supports arrays of constants and variables, and indexed 
by variables (or other expressions). These array expressions are most general. 
However, domain reduction in OPL is weaker than possible, for instance the 
reduction for an index variable depends on its position (PH, p. 100). 

In this work we study constraint propagation enforcing arc-consistency for 
general array expressions. Arrays are multidimensional, and variables can occur 



wherever constants can. An expression x = a[yx, ... ,y n ] is seen as a constraint 
on the variables x, and yi, . . . , j/ n , and all variables collected in the array a. 

Example. Consider an application of arrays. Assume a conventional crossword 
grid, with entries for words in the rows and columns and remaining fields black- 
ened. Further consider a set of words, a subset of which is to be filled into the 
entries in the grid. A natural formulation of this problem as a CSP is to take for 
each word entry a variable Ei whose initial domain is the set of words that fit 
in the entry length- wise. 

The words, split up in their letters, are collected in a two-dimensional array 
/ such that l[w,p] represents the letter of word w at position p. The conditions 
on crossings of entries are then easily stated as constraints. A crossing of field 
Ei at position 4 and field E2 at position 3 is stated as i[£i,4] = Z[i?2,3]. An 
additional alldif f erent constraint on the Ei ensures that no two word entries 
contain the same word. 

Enforcing arc-consistency for array expressions solves some instances of the 
crossword problem without any backtracking ( Jig] , p. 140, which uses a custom 
constraint for crossing entries that is equivalent to the one here). 

2 Preliminaries 

A constraint satisfaction problem (C; T>) consists of a set of variables (implicit 
here), a set T> of domain expressions x £ D x that associate with every variable 
a set of admissible values, and a set C of constraints. A constraint is a relation 
on a set of variables that is a subset of the cartesian product of their domains. 

A solution for a constraint is an assignment of values to its variables that 
is consistent with the constraint. A solution for a CSP is an assignment that 
is a solution for all its constraints. A CSP, or a constraint, is satisfiable if a 
solution exists. A domain value, or a partial solution, is supported if it is part of 
a solution. 

Local consistency notions, weaker approximations of (global) satisfiability, 
are essential in constraint solving. A central one is arc-consistency ([||). We 
disregard the arity of constraints and define: a constraint is arc- consistent, if all 
domain values of all its variables are supported. A CSP is arc-consistent if all 
its constraints are. 

2.1 A Rule-based Formalism 

Constraint programming can be seen as transforming a CSP into one or several 
simpler but equivalent CSPs in a rule-based way. This view allows separate con- 
sideration of the reductive strength of some set of constraint propagation rules, 
and its scheduling. The transformations on CSPs lend themselves to a declara- 
tive formulation. We adopt the proof theoretic formalism of |lj, and introduce 
the elements relevant here. 



A transformation step from a CSP V, the premise, to a CSP Q, the conclu- 
sion, by application of a rule (r) and possibly subject to a side condition (C) on 
V is represented as 



Two CSPs V and Q are equivalent if all variables in V are present in Q and 
every solution for V can be extended to a solution for Q by an assignment to 
variables only in Q. If a rule application preserves equivalence then the rule is 
sound. 

A rule is required to be relevantly applicable, that is, the result Q must be 
different from V in the sense that the set of domain expressions or the set of 
constraints changes. If a rule, or a set of rules, is not applicable to V then V is 
stable or closed under it. Applying a rule to a constraint means applying it to 
the CSP consisting only of this constraint. 



Notation. Domain expressions v € D v used in rules are implicitly represented in 
the set V. Replacing a domain expression present in V is denoted by V, v 6 D™ ew . 
If in V = (C; V) the set of constraints consists of only one constraint, C = {con}, 
then we may just write V = (con;T>). The expression sh^t denotes a substitu- 
tion, assignment, or mapping, from s to t. 



Example. We illustrate these concepts with a rule-based characterisation of arc- 
consistency. A constraint C is arc-consistent if for all variables v of C and all 
values d e D v an instantiation of v to d in C, written C{v^>d}, retains satis- 
fiability. If C{vh^d} is not satisfiable then d is redundant and can be removed 
from D v . The resulting CSP is equivalent to the original one. This principle is 
captured in a rule: 

Lemma 1. A satisfiable constraint C is arc- consistent iff it is closed under the 
application of 

(C • T>) 

(ac) — — — ■ — „ ; : „ ; C{v^d\ has no solution 

□ 



3 Arc-consistency for Array Constraints 

An array a of arity n is a set of mappings index variable, index is a unique 
n-tuple of constants, variable is a variable with a domain. The array expression 
a[bi, . . . , b n ] evaluates to v, if a contains a mapping (bi, . . . , b n )<-^v, otherwise it 
is not defined, (in what follows it is assumed that indices accessing a are valid). 
Note that arrays of constants come as a specialisation of this model. 



3.1 Simple Array Constraints 



Array expressions a[yi, . . . , y n ] are functional. The simplest extension to a con- 
straint is the equality constraint C = (x = a[yi, . .., y n })- We establish arc- 
consistency first for this case, and discuss subsequently compound (nested) array 
expressions. Also, occurrences of variables are restricted in that no variable in the 
constraint may occur more than once (C is a linear). Note that the variables of 
C are x, yi, . . . , y n , and all variables v for a valid (pi, . . . ,b n ) and (61, . . . , b n ) 1— >d 
in a. Such a v will from now on be denoted directly as a[b\, ...,&„]. 

Theorem 2 (Arc-consistency for arrays). A satisfiable linear equality con- 
straint (x = a[y\, . . . , y n }) is arc- consistent iff it is closed under the rule set 
T^arr : 



(arr^) 



(x = a[yx, ...,y n ];T>) 



a[yi,...,y„] ; V,x e D x n {\J bieDy . A>[6i,...,&„]) 



, v {x^a[ yi ,...,y n ];V) 

™ (x = a[y 1 ,...,y n } ] V, yk GD yk \{b}) {vl 



(C y ) : D x nl |J = 

\6« €D Sj , b k =b J 

(arra) (a;^a[yi,...,y n ]; V) 

= a[yi, . . . ,y„] ; X>,a[&i, . . . ,6„] G £> a [bi,. ..,&„] n D x ) 

(C„) : x ... x = {(&!,...,&„)} 

□ 



Proof. 

(<*=) Suppose C = (x = a [yi, . . . , y n ]) is closed under 7\Larr- Then all values in 
domains of variables in C are supported. 

(1) Take some d G D x . C is closed under (arr x ), thus also d G 

(jJb t eD y . Ax[6i,...,6„])- Tncn there exists some (61, ...,&„) with d G A»[&i,...,6„]- 
This index and a [61, . . . , b n ] 1 — >• cZ support x^d. 

(2) For some & G D yfc consider the necessarily failing condition of (arr y ). Thus 
a value d exists in both D x and Aj[6i,...,6 n ]i f° r some (61, . . . , b n ) with &fc = b. 
Assigning the 6, to the j/j, and xi->d and a [61, ... , 6„] h^c£, is a solution supporting 
b. 

(3) Consider a value d G ..,&„]) an d the following cases: 
(3.1) (61,..., b n ) i D V1 x ... x £)„„ 

The index (j/i, . . . , y n ) can not select the variable a[b\, . . . , &„]; however, C re- 
mains satisfiable. Therefore, there is a solution for C that is indifferent to the 
value of a[bi, . . . , &„], and so supports a[&i, . . . , 6„] n->d. 



(3.2) (&!,..., b n ) £ D yi x ... x Z? y „ 

(3.2.1) {(b 1 ,...,b n )} = D yi x...xD yn 

Here the condition of (arr a ) is fulfilled, its consequence holds and with it d £ D x . 
A supporting solution is therefore x>— >d, a[bi, . . . , b n ] i— >d, yi^bi for all i. 

(3.2.2) some D Vk contains more than one element 

Consider some index (b' l: . . . ,b' n ) with b' k ^ bk that also fulfills DxDD^^ b i ] ^ 
0. Such an index exists because otherwise (arr y ) would be applicable. Choose a 
d' £ D x and instantiate md', a[6' 1; . . . , 6^] >— > d', yi^>b[ for all i. This solution 
to C does not assign to a[b\, . . . , 6 n ] and hence supports a\b\, . . . , b n ] t-^d. 

(=>) Suppose here that C is not closed under 7^arr- Then domains of some 
variables in C contain unsupported values. 

(1) Assume (arr x ) is applicable, that is, D x D UbieDy ^a[6i, ...,&„]• Then there is 
some value d £ D x ,d £ O a [6i,...,6„] f° r an (^i> ■ ■ ■ >b n ) £ D Vl x . . . x D Vn . Clearly, 
d is not part of any solution. 

(2) Suppose some bk £ D Vk could be removed by (arr y ). From the condition of 
(arr a ) it follows that with yk<—>bk no index (pi, . . . , b n ) can be found that allows 
the same value for x and a[b\, . . . , b n }. 

(3) For a singleton index domain and so a possible application of (arr a ) consider 
a[b\, . . . , b n ]^d with d D x . Such a d can not be supported by x. □ 

Linearity requirement. It is necessary to restrict occurrences of variables. Con- 
sider the array xor = {(0, 0) i-> 0, {(0, 1) h-> 1, {(1, 0) ^ 1, {(1, 1) h-> 0} and the 
CSP V = (0 = £or[y, y]; {y G {0, 1}}). 7-* is inconsistent but stable under 7?.arr- 

Origin oJTZqxv- Each rule in 7?-arr can be derived as an instance of the general 
rule (ac) in Lemma |l|. Such a derivation, perhaps unsurprisingly, proceeds along 
the same case distinctions as in the (<;=) part of the above proof. We believe the 
derivation to be interesting in its own right, but choose here the proof for its 
relative brevity. 



3.2 Arc-consistency and Compound Expressions 

The following result allows decomposition of nested array expressions and equal- 
ity constraints for the purpose of establishing arc-consistency. Expressions such 
as l[w,p] = l[w' \p'] from the crossword example are decomposed with a fresh 
variable into v = l[w,p] and v = l[w',p'], upon which arc-consistency can be 
enforced independently. 

Lemma 3. Assume Ct = (s = t(v)) and C v = (v = r) be linear constraints on, 
apart from v, distinct sets of variables. The constraint C = (s = t{v<— >r}) is 
arc- consistent if Ct and C v are. □ 

Proof. Suppose Ct and C v are arc-consistent. 

Any solution for Ct assigns a value to v that is also supported by a solution 
to C v , and vice versa. Due to the conditions on variables, such solutions do not 



assign to the same variables. Therefore there union is also a solution for C . Thus, 
a supporting solution for any domain value of a variable in Ct,C v , and C, can 
be extended to a supporting solution for C . 

Hence, C is arc-consistent. □ 



3.3 Domain Reduction and Transformation 

As instances of (ac), the rules in IZ&tt are domain reduction rules by type. From a 
semantical, and particularly from an operational, point of view, however, it may 
be worth to have instead transformation rules which change the representation 
of constraints. 

Consider (arr Q ), which applies if the index is fully instantiated. That means 
we can also dispense entirely with the array look-up: no choice is left. The array 
expression can be replaced by the selected variable. An alternative to (arr a ) 
would thus be 

(x = a[y 1 ,...,y n ];V) 



(x = a[bi, ...,&„]; V,a[bi, ...,b n ] G D [in,...,5 B ] H D x ) 

This rule is now both a transformation rule and a domain reduction rule. Note 
that the domain reduction takes place between variables. In presence of rules 
for primitive equality constraints (a; = y) one can simplify even more into a pure 
transformation rule: 

, M (x = a[y u ...,y n };T>) 
{x = a[bi, . . ,,b n \ ; V) 

The combination of (arrJJ and rules for (x = y) is equivalent to (arr a ). 



4 A Non-naive Algorithm 

An exhaustive application of 7?-arr is computationally expensive, in part unavoid- 
able due to the strength of arc-consistency, and the large number of variables 
involved in array constraints. An inefficiency that can be remedied is the large 
number of set operations on domains, due to fact that individual array variable 
domains D [6i,...,6„] are read and processed multiple times. 

The algorithm arr-ac (Fig. [l]) reads every D a t blt ^ tb i addressable by 
(j/i, . . . , y n ) at most once. Consider T = D x n -D a [i>i,. ..,&„] for some (&i, . . . , b n ) G 
D Vl x ... x D Vn . T is a subset of the intersection in the conclusion of (arr^), 
so it is necessarily part of the new domain of x, and only D X \T instead of D x 
needs to be subjected to further restriction. With regard to (arr v ), a nonempty 
T implies that the side condition fails. Thus, no fefc of (b\, . . . , b n ) can be removed 
from D Vk by (arr y ). 

Note that arr-ac makes a positive guess whether values are supported. If 
indeed in the end some domain is reduced then arr-ac needs to repeat the run. 
Indeed, if at some before the regular end of the run as described in Fig. p] 



for all i: 
for all i: 

X:= D x 
S:=0 



Bi ■= Dy 

Y := D v . 



index 

potentially redundant for yi 
potentially redundant for x 
indices skipped for X 



while B / and some Y k / loop for Y and X 

choose and remove (61, . . . , b n ) from B 
if some bk € Yt then 

T:= D,fl% I,, 

if T / then 

for all i: Y, := ^{6;} 
X := X\T 
else 5:=5U{(6i,...,6„)} 



while 5/0 and X ^ 

choose and remove (61, 
X := X\D a [ 6li ... ii)ri ] 



, b n ) from 5 



for all i: D y 

D x := D x — X 



Yi 



rest loop for X 



remove redundant values 



Fig. 1. arr-ac (core) 



it is definite that some domain D Vj will be reduced, the run could terminate 
immediately, commit the change to D yj , and restart. 

The core part of arr-ac can itself be regarded as a complex domain reduction 
rule, encapsulating (arr^) and (arr^). The rule set {arr-ac:core, (arr a )} establishes 
arc-consistency. 

Example. Consider x G {B, C, D}, yi G {1, 2}, t/2 € {1, 2, 3} and (x = a[y\, 2/2]) 
and let a be defined as the array of constants 



(2/1,2/2) 


1 


2 


3 


1 


A 


B 


C 


2 


D 


E 


F 



The constraint is arc-consistent, which arr-ac verifies as follows. First it reads 
a[l, 1] = A. Nothing is done. It follows a[l, 2] = B. B is in D x , so B is a supported 
value for x, and 1 is supported for y\ and 2 for y^. The next step is reading 
a[l,3] = C. This supports C for x and 3 for y 2 . Finally, a[2, 1] = D is reached. 
This supports the last missing value D for x, and moreover 2 for y\ and 1 for y 2 . 

Support for all values in the domains was found, hence arc-consistency 
is established. Only one incomplete run was necessary, skipping the indices 
(2, 2), (2, 3) that are still permissible by (t/i, y 2 )- □ 

For one run of arr-ac (and ignoring X here), the number of iterations that 
enter the computation of T has an upper bound of O (d n ) with d the maximal size 
of the domains of the yi. This reflects the number of possible different indices 
(61, . . . ,b n ). The lower bound, on the other hand, is only 0{d). It is reached 



when every iteration reduces all (nonempty) Yi by an element, and occurs if the 
constraint is arc-consistent and every instantiation of (j/i, . . . , y n ) is part of a 
solution. 

An operationally useful side effect of arr-ac is that it can also yield the 
variables that contain the supporting values. Initially, all variables a[b\, . . . , b n ] 
are part of the constraint, whereas after complete instantiation of the index 
(j/i, ...,?/„) only the variable a[yx, . . . , y n ] is constrained and contains support, 
arr-ac regards those variables a[b\, . . . , b n ] as supporting, for which the intersec- 
tion T is nonempty. 

arr-ac was implemented in ECL l PS e (Q), using the finite domain primitives 
of lib(fd). An implementation of 7?-arr in the same environment was compared 
to arr-ac by testing it against an instance of the crossword problem, and was 
roughly 50% slower. 

5 Final Remarks 

5.1 Related Work 

The established precursor of array constraints is the element constraint (Q). 
It is the one-dimensional specialisation, and usually the look-up list that links 
index and result is restricted to consisting of constants. 

Arrays in OPL |^]) are similarly general as in this work. In j|] on OPL++ 
a model of the stable marriages problem is described that employs an array of 
variables indexed by a variable. Constraint propagation of array expressions in 
OPL is strictly weaker, however. For all three cases treated by 7^arr we could 
construct simple examples using small 2-dimensional arrays in which reduction 
of domains is possible but not performed, see Figures | and | 

|1 describes an implementation of element using indexicals in AKL(FD), 
in which the look-up list can consist of domain variables. It is equivalent to a 
one-dimensional instance of 7?-arr- 

In j^] a new constraint case is proposed that subsumes multidimensional 
array constraints with arrays of constants. An algorithm which seems similar in 
effect to using 7?.arr, based on graph theory, is outlined. 

P|, on unifying optimisation and constraint satisfaction methods, studies a 
continuous relaxation of element with a look-up list of variables with continuous 
domains by using a cutting-planes approach. 

5.2 Conclusions 

We study here the use of arrays in constraint programming mainly from a theo- 
retical point of view. There are good arguments suggesting that arrays are bene- 
ficial in constraint models, however. Indices on objects are basic in mathematics, 
element is implemented in many constraint systems. Arrays with multiple di- 
mensions have long been used in imperative, now object-oriented, languages. 
These language styles obviously inspired OPL (|ll|, more so ||), a successful 



constraint programming system. Yet it would be desirable to have large exam- 
ples of uses of multidimensional arrays. 

Such problems could also be used to evaluate the use of arc-consistency as the 
objective in constraint propagation. It is now clear from practical experience that 
the notion of consistency that is most advantageous depends on the problem. 
Sometimes a weaker notion such as bound or range consistency might suffice, 
for example applied in the early stages of solving and later replaced by full arc- 
consistency. 7^arr provides a starting point for obtaining reduction rules for those 
consistency notions, which are subsumed by, yet very similar to, arc-consistency. 
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// 
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1, m } 


{ 1 } 



Fig. 2. OPL: non-applied (arr x ), (arr^) 



enum Dy { i , j , k } ; 

enum Da { p, q, r >; 

var Da a [Dy] ; 

var Da x ; 

var Dy y; 

solve { y = j ; 

x <> q; // OPL arc-consistency 

x = a[y] ; // a[j] in { p, q, r } { p, r } 

}; 



Fig. 3. OPL: non-applied (arr a ) 



